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ANALYSIS OF PLASMAS GENERATED BY FISSION FRAGMENTS 


Jerry E. Deese and H. A. Hassan 
North Carolina State University 

SUMMARY 

A kinetic model is developed for a plasma generated by fission fragments 
and the results are employed to study He plasma generated in a tube coated with 
fissionable material. Because both the heavy particles - and electrons play im- 
portant roles in creating the plasma, their effects are considered simultaneous- 
ly. The calculations are carried out for a range of neutron fluxes and pres- 
sures. In general, the predictions of the theory are in good agreement with 
available intensity measurements. Moreover, the theory predicts the experi- 
mentally measured inversions. However, the calculated gain coefficients are 
such that lasing is not expected to take place in a helium plasma generated by 
fission fragments. The effects of an externally applied electric field are 
also considered. 


INTRODUCTION 

Increased interest in gas core reactors’* - and the recent demonstration of 
2-5 

direct nuclear pumping focused attention on plasmas generated by the high 
energy fission fragments. Such systems are rather complex and the plasma gen- 
erated in them is, in general, not in thermal equilibrium. Therefore, before 
one can predict their performance characteristics, one needs to develop a de- 
tailed and self-consistent kinetic model capable of predicting the behavior of 
the plasmas generated in these devices. 



Several authors have analyzed the space-dependent volumetric production of 
ions by fission fragments passing through a background gas. Both Leffert, 
et al.^ and Nguyen and Grossman^ derived expressions for the spatial distribu- 
tion of fission fragment ion production. In both of these analyses an energy 
independent empirical value W is assumed for the amount of energy required to 
produce an ion pair. The former utilized both linear and square energy-loss 
models for the heavy particles, while Nguyen and Grossman used the Bohr stop- 
ping equation for fission fragments. Rather than rely on an empirical constant 
8 9 

Miley and Thiess ’ derived an expression for the ionization and excitation 
rates which takes into account the effect of the energy distribution of the in- 
cident particles. Such a calculation requires estimation of the energy depen- 
dent cross sections for excitation and ionization by heavy charged particles. 

A Bethe-Born type representation was employed for the case of helium excitation 

by alpha particles and fission fragments^ and later Guyot, et al.^ employed 

11 

Gryzinski cross sections for helium ionization by alpha particles and lithium 
ions. 

The calculation of the electron energy distribution function in electric 

discharges and in the absence of a high energy volumetric source is a standard 

procedure. However, there are only a few analyses of distributions resulting 

from a flux of high energy particles where there is a high energy volumetric 

electron source. Calculations using a Monte Carlo method were carried out by 
12 13 

Wang and Miley . Later, Lo and Miley used a simplified version of the 
Boltzmann equation to determine the electron energy distribution in a helium 
plasma produced by a mono -energetic electron source. More recently, Hassan and 
Deese"^ presented a more elaborate Boltzmann equation formulation which took 
into consideration the primary electron spectrum. 


2 



In general, theoretical studies of excited state densities have assumed 
the electron distribution function to be a Maxwellian at some characteristic 
temperature. Russell used such an assumption in calculating excited state 
densities in argon. Leffert, et al.^ and Rees, et al."^ also studied noble gas 
plasmas assuming the electron distribution function to be Maxwellian. More re- 
cently Maceda and Miley^ calculated the number densities of the helium excited 

13 

states using the non-Maxwellian distributions of Lo and Miley ; their results 

indicate a number of possible inversions. 

In addition to the above analyses, there exists a number of experimental 

investigations dealing with nuclear pumped lasers and laser enhancement. Work 

carried out before 1972 is summarized in reference 18 while a summary of the 

nuclear laser effort at the University of Illinois along with an exhaustive 

list of references on virtually every aspect of radiation produced noble gas 

19 

plasmas is included in the work of Thiess . Of particular interest here are 
experiments studying individual atomic transitions at various pressures and ad- 
ditive concentrations under fission fragment excitation. The earliest study of 

20 

fission fragment excited spectra is that of Morse, et al. who examined the 

21 

effects of fission fragment radiation on He, Ar, N^, and air. Guyot measured 

the production of helium metastables by B^(n,a) fission fragments, while 
22 

Walters measured the relative intensities of the various transitions in both 
helium and argon. 

As a result of the numerous research efforts outlined above, actual nuclear 

pumped lasing has just recently been reported by several authors. McArthur and 
2 

Tollefsrud reported lasing action in carbon monoxide as a result of nuclear ex- 

3 

citation only. Helmick, et al. demonstrated direct nuclear pumping in He-Xe 

4 

gas mixtures. A third case of direct nuclear pumped lasing is that of DeYoung 


3 



in a neon-nitrogen mixture. All of the above lasers employed tubes coated with 

fissionable material. The work of Jalufka, et al. ^ employed a volume source of 

3 3 3 

fission fragments; the He(n,P) H reaction was employed to excite a He-Ar 

laser. Obviously the nuclear pumped gas laser research effort is still in its 
early stages with the validity of the concept having been demonstrated only re- 
cently. 

The object of this investigation is the development of a theoretical model 
for a plasma generated by high energy fission fragments. The kinetic model 
treats particles in different quantum states as different species and uses the 
multifluid conservation equations of mass, momentum and energy to describe the 
resulting system. It takes into consideration the following kinetic processes: 
ionization, excitation and deexcitation, radiative recombination, spontaneous 
emission, associative ionization and dissociative and collisional recombination. 
Because both the heavy particles and electrons play important roles in creating 
the plasma, their effects are considered simultaneously. The rates of reactions 
involving electrons were calculated using electron distribution functions ob- 
tained from a solution of a Boltzmann equation appropriate for plasmas generated 
by fission fragments'^. 

The above model is employed to study a helium plasma generated by fission 

fragments. Helium was chosen because of the availability of experimentally 

22 

measured cross sections and rates and in-reactor measurements . In general, 
the results show good agreement with experiment. Moreover, they indicate a 
number of possible laser transitions; all of them, however, are in the IR re- 
gion. 
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ANALYTICAL FORMULATION 


The systems to be modeled here are those appropriate for nuclear pumped 
lasers. Typically, they consist of tubes coated with fissionable material and 
filled with gas at some given pressure and temperature. The tube is placed in 
the high neutron flux region of a reactor. Under neutron bombardment fission 
fragments emerge from the coating and enter the gas. The ensuing energy trans- 
fer results in ionization and excitation of the background gas. A schematic of 
the slab geometry employed in this analysis is shown in figure 1. 

Treating particles in different excited states as different species one 
can utilize the multifluid conservation equations to describe the plasma gen- 
erated by the fission fragments. For the conditions under consideration, the 
steady state approximation is appropriate. In this approximation, the effects 
of gradients are assumed negligible. Thus, the conservation of species equa- 
tions. 


8n -* 

-x- 5 - + V * (n U ) = I + R 
dt s s s s 


reduces, as a result of this approximation, to 


(i) 


I + R = 0 . (2) 

s s v 

In the above equations, s is a charged particle or any quantum state of the 

background gas and, for species s, n is the number density, U is the velocity 

s s 

and I g and R g are the production rates per unit volume resulting from nuclear 
and kinetic sources, respectively. 

For very low Mach numbers the momentum equation reduces to 

p - const. (3) 

where p is the pressure, while the energy equation takes the form 

p ± = Q <A) 
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where P_^ is the power input and Q is the conduction and radiation losses. For 
an optically thin highly conducting gas, equation (4) reduces to 

T - T ± (5) 


where is the initial gas temperature. 

The properties of the electrons are determined from an electron Boltzmann 

equation. The equation employed is that developed in reference 14. For the 

high pressures of interest the plasmas generated by the fission fragments are 

slightly ionized. Therefore, using the Lorentz gas approximation, the result- 

14 

ing Boltzmann equation for a quasi-steady plasma can be written as 


v 9 
3 3x . 


9f 
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v 9x. 
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mv 9v 
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- y [v' 2 Q (v')f ' - V 2 Q (v)f ] + Of /at) = 0 

Lt * n R n n 


( 6 ) 


where 


1 ,2 1 2 . 1 2 
-y mv = t mv + r mv 
2 2 2 s 


(7) 


v is the velocity, e is the electric charge, m is the electronic mass, M is the 

mass of the heavy particles, N is the gas number density, V Q is the collision 
2 

frequency, 1/2 mv g is the excitation energy, Q g is the excitation cross sec- 
tion and (9f /9t) is the source term resulting from primary and secondary ion- 

o c 

ization and recombination. An explicit expression for (9f /9t) together with 

° c 

the method of solution of equation (6) is given in reference 14. 


The quantities I and R g must be determined before the above system of 
equations can be solved for the various species present. To determine I and 
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R g , one needs to specify the important kinetic processes in the system. The 
major reactions included here, which are appropriate for noble gases, are 

1. Fission fragment excitation 

ff + X -> X. + ff 
3 

2. Fission fragment ionization 

ff + X -»• X + + e + ff 

3. Spontaneous emission 

X. X. + hv 
i 3 

4. Ionization and recombination 

X i +eJX + +e+e 

5. Electron excitation and de-excitation 

X. + e t X. + e 
i 3 

6. Radiative recombination 

X + + e X + hv 

7. Excitation transfer 

X + X t Xj^. + X 

8. Dissociative recombination 

X* + e -> X ± + X 

9. Collisional recombination 

X + + X + X + X2 + X 

10. Associative ionization 

X + X ± -► X* + e 
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11. Electron stabilized recombination 

X 2 + e + e->X i + X+ e 

12. Neutral stabilized recombination 

xt+e+X->X. + X+ X 
2 i 

13. Metastable ionization 

X* + X* X + + X + e 

14. Charge transfer 

X + + M -*■ M + + X 

15. Penning ionization 

X* + M -> M + + X + e ( 8 ) 


where M represents a substance other than the background gas, or an impurity. 

The term I is a result of reactions of Type 1 and 2 while R g is obtained 
from reactions of Type 3-15. In general, for a reaction of the type 


a. A. + a. A. •> a A +a^A^ 
xi 33 ss tt 

the contribution to the production rate of species s is 


(9) 


a K 


a. a. 


( 10 ) 


where a^(Jl = i, j, s, t) denotes a stochiometric coefficient. The quantity K 


is the forward rate coefficient and is given by 

K = | f . f . g. . a. . dV. dV. 

l j 13 l J 


( 11 ) 


where, for species i, f^ is the energy distribution function, is the veloc- 

• ->■ i 

ity* g_jj = | V± — I and CTy is the collision cross section. If i represents 
the stationary background gas and j a fission fragment or an electron, then 


equation ( 11 ) reduces to 


8 


( 12 ) 


k = f. a v. dv. . 

J J 3 3 

The rate coefficients for reactions involving the background gas (or gases) 
are usually obtained from experiment. On the other hand reactions involving 
the fission fragments and electrons can, in principle, be calculated according 
to equation (12) from collision cross sections determined from experiment or 
theory and appropriate distribution functions. In this work the electron dis- 
tribution function is calculated from equation (6) . Unfortunately, the situa- 
tion with regards to fission fragments is not well understood; this is because 
the fission fragments are characterized by initial energies ranging from 50 to 
115 Mev, initial charges from 16e to 24e and masses from 70 to 160 atomic mass 
units. In addition no data is available on cross sections for ionization and 
excitation by fission fragments. Because of these uncertainties, the contribu- 
tion of the fission fragments to excitation and ionization was estimated using 
two different approximate methods. In the first, the procedure outlined in re- 
ference 6 was used to estimate the average energy deposited in the gas per unit 
volume per unit time, E^. The average number of ions produced per unit volume 
per unit time is given by dividing by W, the energy expended per ion pair 
produced. Similarly, the total number of excited states produced is determined 
by dividing E,. by W , the energy expended per excited state produced. Rees, 

IT GX 

16 

et al. determined that the total excited particle production rate from fis- 
sion fragments is .53 times the ion production rate, thus 

W = W/.53 . (13) 

ex 

This procedure determines only the total number of excited states produced by 
fission fragments and some model for the distribution of these states must be 
adopted. Because of the absence of a generally accepted procedure for the 



distribution of excited states a number of models have been employed here and 
these are discussed under Results and Discussions. 

9 

The other approach uses the procedure of Thiess and Miley , or Guyot, et 
al.^, which is based on a heavy particle distribution function derived from a 
semiempirical slowing law together with ionization and excitation cross sec- 
tions based on the Born or the Gryzinski approximations. To utilize this pro- 
cedure one needs to assume that the fission fragments fall into two groups: a 

light group with an average mass number of 96, an average charge of 20e and an 
average initial energy of 98 Mev; and a heavy group with an average mass number 
of 140, an average charge of 22e and an average initial energy of 67 Mev. A 
convenient summary of all formulas needed to calculate the contribution of the 
fission fragments to ionization and excitation by the two methods outlined 
above is given in reference 10. 


RESULTS AND DISCUSSION 

The above model is applied to a study of He plasma generated by fission 

22 

fragments and the results are compared with the measurements of Walters who 
employed a tube of radius 1.85 cm coated with UgOg. To conform to the condi- 
tions of the experiment, the calculations allow for the presence of a nitrogen 
impurity in the system. For a given pressure, temperature, neutron flux and 
tube dimensions, the above model is capable of predicting the number densities 
of the helium excited states, the atomic and molecular ions and the electrons. 
From this, one can calculate the relative intensities and gain coefficients. 

For the calculations presented here, the gas temperature is assumed constant at 
300°K, while the pressure ranges from 100-760 Torr (1 Torr = 133 N/m^) , the 


10 


neutron flux from 3.8 x 10^ - lO 1 ^ neutrons /cm?sec . and the nitrogen con- 
centration from .001 to 50 parts per million. 

All helium excited states with a principal quantum number of 5 or less are 
included in the calculations. The rates or cross sections for the reactions of 
Types 3 through 15 indicated in equation (8) were obtained from references 23- 
54 while Einstein coefficients for spontaneous emission were taken from refer- 
ences 55 and 56. A listing of the reactions included in the analysis with re- 
ferences to cross section and rate data is given in Appendix A. Because the 
experimental data is incomplete the products or the rates of some reactions had 
to be estimated and these are discussed next . 

At the high pressures of interest here the recombination process in noble 
gases is complicated by the formation of molecular ions. For pressures greater 
than 5 Torr the reaction 


He + 2 He -*■ He* + He 

23 


(14) 


quickly converts atomic ions into molecular ions . The molecular ion recom- 
bination is governed by three reactions^ ’ ^ 


e + He 2 ->■ He + He 

e + e + He* + He + He + e 

e + He + He* + 2 He + He (15) 

* 

where He denotes an excited state. The distribution of excited states produced 

25-27 

by these reactions is not well known. However, recent studies indicate 

that at least 70% of the excited states produced are atomic metastables in the 

pressure range of interest here. Although potential energy curves indicate 

3 

other states such as the 2 P states are produced in molecular ion 


11 


recombination' 


57 * 

, it is assumed that He in the reactions indicated in equation 
(15) are atomic metastables. 

Rates for associative ionization are available for atomic states of prin- 
cipal quantum number 3 and triplet states of principal quantum number 4, refer- 
ences 26, 28-30. Cross sections for the excitation transfer reactions are 
available for n = 2 and n = 3 but only for the triplet states when n = 4, re- 
ference 26. Because these processes are important in the determination of the 
final distribution of excited states some estimate of the unavailable rates is 
required. It is assumed here that the associative ionization rates for n = 5 
states are the same as the corresponding n = 4 states while the associative 
ionization cross sections for the n 1 F states are assumed equal to those for the 
n^D states. Moreover, three calculations were carried out in which the cross 
sections for associative ionization and excitation transfer reactions of the 
triplet n = 4 states were assumed to be one-third of, equal to, and three times 
the corresponding singlet reaction cross sections. Comparison of calculated 
and measured excited states showed that better agreement with experiment is ob- 
tained when Q(4 L X) < Q(4 3 X). Therefore, unless indicated otherwise, all calcu- 
lations reported here employ this assumption. 

As indicated earlier, when a 'W' value is used to estimate the total num- 
ber of excited states further assumptions are needed to indicate which states 

are excited. From a study of the spectra of noble gases generated by the im- 

58 

pact of alpha particles, Bennett suggested that most of the excited states in 

He are nT* states with the majority of them in the 2 1 P state. His conclusion 

was based on the fact that calculated excitation cross sections, based on the 

Born approximation, are highest for the n^P states. This assumption is con- 

21 

trasted with that of Guyot where he assumed that most of the excited atoms 
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resulting from the impact of alpha particles are in the 2^S state and these are 
3 

converted to the 2 S states by electron impact. Calculations using both of the 

above models were carried out; for these calculations a value of W =85 ev 

’ ex 

was employed. 

22 

Comparison of the results with the experiments of Walters requires that 
the effect of a nitrogen contaminant be taken into consideration. The gas used 
in the experiment contained a nitrogen contaminant of approximately 50 parts 
per million. The effect of the pressure on the number densities of the elec- 
trons and the two metastable states is given in figure 2 while figure 3 shows 
the effect of pressure on several of the higher states with and without a 50 

parts per million nitrogen contaminant. For these figures, the neutron flux is 
11 2 

3.8 x 10 neutrons /cm. sec. Considering pure helium first, it is seen that the 

electron number density decreases as the pressure increases. This is because 
the dominant helium recombination process is the neutral stabilized recombina- 
tion of the molecular ion and this process becomes more efficient as the pres- 
sure increases. The dominant processes governing the excited states number 

densities are electron excitation from the ground state and the metastable 
3 

state 2 S, excitation transfer, associative ionization and spontaneous emission. 
As the pressure increases, the decrease in the electron density coupled with 
the increased effectiveness of the associative ionization and excitation trans- 
fer results in the decrease of the excited states shown in figure 3. Only the 

2"^S and 2 L P states show an increase as the pressure increases; the rapid recom- 

3 

bination resulting from the last reaction in equation (15) accounts for the 2 S 
behavior . 

As is seen from figures 2 and 3, the effect of nitrogen on the various spe- 
cies is quite significant . Both the helium molecular ion and the metastables 
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are quickly converted to nitrogen ions whose dominant recombination process is 
dissociative recombination. This not only lowers the concentrations of the He 
ions and metastables but changes the electron number density variation with 
pressure as well. When nitrogen is present the dominant recombination process 
is the two-body dissociative recombination. This process does not become more 
efficient as pressure increases and as a result the electron number density in- 
creases with an increase in pressure. Thus, the net effect of nitrogen is two- 
fold: it lowers the number density of electrons and metastables making elec- 
tron excitation from metastable states less important and, it changes the num- 
ber density dependence on pressure. As a result of this, one would expect a 
significant effect on the higher excited states as well. Because the electron 
number density increases with pressure the electron excitation rates increase 
with pressure. Thus, as is seen in figure 3, the marked decrease in excited 
state population with pressure in the case of pure helium is not as pronounced 
when nitrogen is present. The figure also shows that the magnitudes of the 
number densities decrease with pressure. 

The results indicated in figure 3 do not exhibit the peak around 200 Torr 

22 19 

shown in figure 4-9 of reference 17. Both Walters and Thiess suggest that 
the associative ionization process may be a three-body process of the form 

He + 2 He -> He* + He (16) 

rather than the two-body process 

He + He ^ He* + e (17) 

assumed here. If this process is indeed three-body in nature its effectiveness 
will increase with pressure. 
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A study of the variation of the number densities with neutron flux for the 
case of a pure helium gas at 100 Torr has also been carried out. The behavior 
of the electrons and metastables is given in figure 4 and that of some higher 
states is shown in figure 5. The number densities of the higher states in- 
crease linearly with the flux at lower neutron flux levels. As the electron 
number density increases, the electron stabilized three-body recombination be- 
comes important. Thus, population densities highly dependent on electron ex- 
citation rates do not increase as rapidly at higher flux levels. 

Because the nitrogen contaminant has a significant effect on the number 
densities, a determination, of the concentration range where its effect becomes 
insignificant is of interest. Figure 6 shows a plot of the electron and meta- 
stable number densities vs. nitrogen concentration at 100 Torr and a neutron 
11 2 

flux of 3.8 x 10 neutrons/cm. sec. As is seen from the figure, the effect of 

“8 

nitrogen is negligible below a concentration of 10 

All of the above results assume the 2^P state to be the only excited state 
produced by fission fragments. Calculations using the procedure of reference 
9, the TM model, have also been carried out. Boltzmann plots showing log 
Xl/g^A = log[hcn^/4Trg^] , where X is the wave length, I is the intensity, A is 
Einstein's transition probability, g_^ is the degeneracy, h is Planck's constant 
and c is the speed of light, vs. £ u , the energy of the upper level, for pres- 
sures of 100 and 760 Torr are given in figures 7 and 8 for the 2 1 P model and in 
figures 9 and 10 for the TM model. Error bars in these graphs indicate the 
spread in Walters' experimental data. Using the TM model, the triplet states 
are considerably underpopulated relative to the singlet states. Because 
Walters found the triplet states to be of the same order of magnitude for states 
with n = 3, it appears that the 2^P model is more appropriate. 
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As is seen from figure 7, the 2 P model gives good agreement with experi- 

3 

ment at 100 Torr; the only discrepancy being the n D states which are overpopu- 
lated compared with the experimental measurements. At 760 Torr the agreement 
of the 2^P model with experiment is not as good as at 100 Torr. All triplet 

3 

states except the 3 S are predicted to have a population greater than experi- 
mentally measured values. This could be due to an improper pressure dependence 
for associative ionization as discussed earlier. In general, singlet states 
have much higher emission coefficients and as a result, cascade losses play a 
larger role in the singlet system. The triplet system on the other hand is 
predominantly governed by associative ionization and excitation transfer pro- 
cesses and these have a greater dependence on pressure. Thus, it is expected 
that triplet states number densities will decrease more rapidly with pressure. 
Calculations have also been carried out assuming fission fragments produce only 

2'*"S states. However, for this case the 2^S number density is higher than the 

3 21 

2 S number density and this is not in agreement with available experiments 

When studying pure helium, two population inversions have been found 

throughout the entire range of pressures, neutron fluxes and fission fragment 

11 11 

excitation models examined on the 3 P - 3 D (95.76y) and 4 P - 4 D (216y) 

22 59 

lines ’ . In addition, for pressures less than 200 Torr and neutron fluxes 

14 2 1 

less than 10 neutrons/cm. sec. the TM model predicts an inversion of the 3 S - 

1 ° 

2 P line (7281 A) . The expression for the gain coefficient is given in Appen- 
dix B. The gain coefficients together with other inversions, operating condi- 
tions and fission fragment models employed are summarized in Table 1. The ad- 
dition of the nitrogen did not change the above results. However, the gain co- 
efficients are slightly decreased because of greater depopulation of excited 
levels . 
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TABLE 1 GAIN COEFFICIENTS 


Transition 

Wavelength 

Excitation 

Model 

Pressure 

Torr 

Neutron 

Flux 

”i. 2 /n r 

Gain 

Coefficient 

3 1 ? - 34 

95.76 ]X 

24 

100 

3.8 x 10 11 

0.0 

0.28 x 10" 9 

34 - 3 1 D 

95.76 y 

24 

760 

3.8 x 10 11 

0.0 

0.76 x 10" 11 

3 X P - 34 

95.76 y 

24 

100 

1.0 x 10 16 

0.0 

0.14 x 10 -5 

3 1 ? - 34 

95.76 y 

24 

760 

3.8 x 10 11 

5 x 10 -5 

0.42 x 10' 11 

3 1 ? - 3-4) 

95.76 y 

TM 

760 

3.8 x 10 11 

0.0 

0.15 x 10~ 9 

4 1 ? - 

216 y 

24 

100 

3.8 x 10 11 

0.0 

0.12 x 10 -9 

44 - 44 

216 y 

24 

760 

3.8 x 10 11 

0.0 

0.24 x 10" 11 

4 X P - 44 

216 y 

24 

100 

1.0 x 10 16 

0.0 

0.16 x 10' 6 

4 X P - 44 

216 y 

24 

760 

3.8 x 10 11 ' 

5 x 10~ 5 

0.14 x 10" 11 

44 - 44 

216 y 

TM 

760 

3.8 x 10 11 

0.0 

0.11 x 10 _1 ° 

4 3 D - 4 3 P 

43.94 y 

24 

100 

1.0 x 10 16 

0.0 

0.20 x 10 -5 

3 1 S - 24 

7281 A 

TM 

100 

3.8 x 10 11 

0.0 

0.24 x 10“ 8 

















Figures 7-10 show only states of principal quantum number four or less; 
these states were the only states considered in the determination of the inver- 
sions present. Excitation transfer cross sections are not available for the 
n = 5 states and as a result these states are overpopulated. The importance of 
these reactions can be seen by performing calculations where these reactions 
are neglected for all states. The Boltzmann plot for such a case is shown in 
figure 11. Comparing figures 7 and 11, it is clear that excitation transfer 
plays an important role in determining the relative population of excited 
states and consequently, the presence of inversion. Similar conclusions hold 
when one assumes 


Q(4 1 X) = 3 Q(4 3 X) , 


(18) 


as is seen from comparing figures 7 and 12. 

Calculations have also been carried out in the presence of an externally 

applied electric field of 10 V/cm and the results are shown in figure 13. For 

this calculation the pressure is 100 Torr, the neutron flux is 3.8 x 10 3 ' 3 ' neu- 
2 

trons/cm.sec. and the background gas is pure helium. As is seen from the figure 

the n^P - n^D inversions present in the absence of an electric field disappear 

3 3 

and a new inversion 4 D - 4 P (43.94y) appears. The gain coefficient for this 

-2 

transition is 7.5 x 10 which reflects the large increase in the transfer of 

14 

electron energy to the excited states . It should be noted that this result 
is for the particular conditions specified above and no attempt has been made 
to study systematically the effect of an externally applied electric field on a 
nuclear pumped laser. 

Recently, Jalufka, et al. 3 demonstrated a volumetric nuclear-pumped laser 
3, 3 3 

using TIe(n,P) H reaction to excite a He - Ar laser. The method presented 
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here can, with slight modification be used to analyze this new mode of nuclear 
pumping . 


CONCLUDING REMARKS 

Considering the incompleteness of available cross sections, the predictions 
of the model are in good agreement with available relative intensity measure- 
ments and measured inversions. However, more complete rate data are needed for 
the accurate predictions of other possible inversions involving the higher 
states. The calculated gain coefficients are so small that, because of cavity 
losses, lasing is not expected to occur in helium. Thus, helium is not a good 
candidate for a nuclear pumped laser. The results also suggest the desirabil- 
ity of studying the effects of an externally applied electric field on nuclear 
pumped lasers. 
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APPENDIX A 


REACTIONS INCLUDED IN THE ANALYSIS 

A listing of the reactions included in the analysis with references to 
cross sections and rate data is given below. If a reference is not indicated 
then the reaction is treated in the manner outlined under Results and Discus- 
sion. The principles of microscopic reversibility and detailed balancing are 
used, where appropriate, to relate forward and backward cross sections and 
rates, respectively. 

REACTION REFERENCES 

Type 3: X. X. + hv 55,56 

13 

Type 4: X^ + e £ X + + e + e 


(1) He + e t He + e + e 31 

Type 5 : X^ + e J X^ + e 

(2) He + e J HeC2 1 S) + e 32,33,34,35 

(3) He+e? He(3 1 S) + e 36 

(4) He + e t He'(4 1 S) + e 36,37 

(5) He + e t He(5 1 S) + e 36,37 

(6) He + e % He(2 1 P) + e 38 

(7) He + e £ HeO 1 ?) + e 36,37 

(8) He + et He(4 1 P) + e 36,37 

(9) He + e t HeO 1 ?) + e 37,39 

(10) He + e i He(3 1 D) + e 36,37 
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(11) 

He + e £ He(4^T)) + e 

36,37 

(12) 

He + e He(5 1 D) + e 

36,37 

(13) 

He + e t He(4 1 F) + e 

40 

(14) 

He + e He(5^F) + e 

40 

(15) 

He + e ^ He(2 3 S) + e 

32,33,34 

(16) 

He + e t He(3 3 S) + e 

36 

(17) 

He + e t He(4 3 S) + e 

36,37 

(18) 

He + e £ He(5 3 S) + e 

36 

(19) 

He + e 2 He(2 3 P) + e 

38 

(20) 

He + e t He (3 3 P) + e 

36,37 

(21) 

He + e £ He(4 3 P) + e 

36 

(22) 

He + e X He(5 3 P) + e 

39 

(23) 

He + e t He(3 3 D) + e 

36,37 

(24) 

He + e ^ He(4 3 D) + e 

36,37 

(25) 

He + eJ He(5 3 D) + e 

36 

(26) 

He + e £ He(4 3 F) + e 

40 

(27) 

He + e t He(5 3 F) + e 

40 

(28) 

He (2 3 S) + e £ He(2 1 S) + e 

41 

(29) 

He(2 3 S) + e £ He(3 1 S) + e 

40 

(30) 

He(2 3 S) + e Z He(4 1 S) + e 

40 

(31) 

He(2 3 S) + e. Z He(5 1 S) + e 

40 
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(32) He(2 3 S) + e£ He (2^) + e 40 

(33) He(2 3 S) + e £ HeO 3 ?) + e 40 

(34) He(2 3 S) + e ;£ He(4 1 P) + e 40 

(35) He(2 3 S) + e 2 He(5 1 P) + e 40 

(36) He (2 3 S) + e 2 HeoV) + e 40 

(37) He(2 3 S) + e £ He(4 1 D) + e 40 

(38) He(2 3 S) + e £ He(5 1 D) + e 40 

(39) He(2 3 S) + e £ He(4 1 F) + e 40 

(40) He(2 3 S) + e£ He(5 1 F) + e 40 

(41) He(2 3 S). + e £ He(3 3 S) + e 40 

(42) He(2 3 S) + e £ He(4 3 S) + e 40 

(43) He (2 3 S) + e 2 He(5 3 S) + e 40 

(44) He(2 3 S) + e £ He(2 3 P) + e 42 

(45) He(2 3 S) + e % He(3 3 P) + e 40 

(46) He(2 3 S) + ej He(4 3 P) + e 40 

(47) He(2 3 S) + e t He(5 3 P) + e 40 

(48) He(2 3 S) + ej He(3 3 D) + e 40 

(49) He(2 3 S) + ej He(4 3 D) + e 40 

(50) He(2 3 S) + e t He(5 3 D) + e 40 

(51) He(2 3 S) + e t He(4 3 F) + e 40 

(52) He(2 3 S) + e $ He(5 3 F) + e 40 
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(53) 

He(2 3 P) + e£ He(2 1 P) + e 

42 

(54) 

He(2 1 S) + ei He(2 1 P) + e 

42 

(55) 

He(3 1 S) + e t HeO 3 ?) + e 

30,43 

(56) 

He(3 1 S) + et HeO 3 ?) + e 

30,43 

(57) 

He(3 3 S) + e£ He(3 3 P) + e 

30,43 

(58) 

He(3 3 S) + e t He(3 3 D) + e 

30,43 

(59) 

He(3 1 P) + e| He(3 1 D) + e 

30,43 

(60) 

He(3 3 P) + e t He(3 3 D) + e 

Type 6 : X + + e -> X + hv 

30,43 

(61) 

He + e -*■ He + 24.586 ev 

Type 7: X. +X^ 5^ + X 

44,45 

(62) 

He(3 1 P) + He t He(3 X S) + He 

30,43 

(63) 

He (3 3 P) + He t He(3 3 S) 4- He 

26,30,43 

(64) 

HeO 1 ?) + He Z He(3 1 D) + He 

30,43 

(65) 

He(3 3 P) + He t He(3 3 D) + He 

26,30,43 

(66) 

He(3 3 S) + He Z He(2 3 P) + He 

26 

(67) 

He(3 1 S) + He t He(2 1 P) + He 

26 

(68) 

He (4 3 S) + He Z He(3 3 S) + He 

26 

(69) 

He(4 3 S) + He Z He(3 3 P) + He 

26 

(70) 

He(4 3 S) + He Z He(3 3 D) + He 

26 

(71) 

He(4 3 P) + He t He(3 3 S) + He 

26 
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(72) 

He(4 3 P) 

+ 

He 

(73) 

He(4 3 P) 

+ 

He 

(74) 

He(4 3 D) 

+ 

He 

(75) 

He(4 3 D) 

+ 

He 

(76) 

He(4 3 D) 

+ 

He 

(77) 

He(4 3 F) 

+ 

He 

(78) 

He(4 3 F) 

+ 

He 

(79) 

He(4 3 F) 

+ 

He 

(80) 

He(4 1 P) 

+ 

He 

(81) 

He (5^) 

+ 

He 

(82) 

He(2 1 S) 

+ 

He 

(83) 

He(4 1 S) 

+ 

He 

(84) 

He(4 1 S) 

+ 

He 

(85) 

He(4 1 S) 

+ 

He 

(86) 

He(4 1 P) 

+ 

He 

(87) 

He(4 1 P) 

+ 

He 

(88) 

He(4 1 P) 

+ 

He 

(89) 

He(4 X D) 

+ 

He 

(90) 

He(4 1 D) 

+ 

He 

(91) 

He(4 L D) 

+ 

He 

(92) 

He(4 1 F) 

+ 

He 


X He(3 3 P) + He 
i He(3 3 D) + He 
X He(3 3 S) + He 
X He(3 3 P) + He 
X He(3 3 D) + He 
X He(3 3 S) + He 
X He(3 3 P) + He 
X He(3 3 D) + He 
X He(4 1,3 F) + He 
X He(5 1,3 F) + He 
X He(l 1 S) + He 
X He(3 1 S) + He 
X HeO 1 ?) + He 
X HeCS 3 !)) + He 
X HeO^) + He 
X HeO 1 ?) + He 
X HeO^D) + He 
X He(3 1 S) + He 
X HeO 1 ?) + He 
X HeO^D) + He 
X He(3 1 S) + He 


26 

26 

26 

26 

26 


46,47 

46,47 

48 
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1 


(93) He(4^"F) + He t HeO 1 ?) + He 

(94) He(4 1 F) + He t HeO^) + He 

Type 8: X* + e Z X ± + X 

(95) He* + e -* He(2 3 S) + He 

(96) He^ + e ->- He(2 1 S) + He 

Type 9: X + + X + X^X 2 + X 

(97) He + 2 He He* + He 

Type 10: X + X ■+■ X* + e 

(98) He(3 1 S) + He -> He* + e 

(99) He(3 1 P) + He -*■ He* + e 

(100) He(3 1 D) + He + Re* + e 

(101) He(3 3 S) + He He* + e 

(102) He(3 3 P) + He -> He* + e 

(103) He (3 3 D) + He -> He* + e 

(104) He(4 1 S) + He -* He* + e 

(105) He(4 1 P) + He + He* + e 

(106) He(4 1 D) + He -> He* + e 

(107) He(4 1 F) + He -*• He* + e 


25,26,27 

25,26,27 

23 

28,29,30 

28,29,30 

28,29,30 

26,28,29,30 

26,28,29,30 

26,28,29,30 
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I 


(108) He(4 3 S) + He He* + e 

(109) He(4 3 P) + He + HeJ + e 

(110) He(4 3 D) + He -*■ He* + e 

(111) He(4 3 F) + He -»■ He* + e 

(112) He(5 1 S) + He + He* + e 

(113) He(5 1 P) + He He* + e 

(114) He(5 1 D) + He + He* + e 

(115) He(5 1 F) + He -> He* + e 

(116) He(5 3 S) + He He* + e 

(117) He(5 3 P) + He > He^ + e 

(118) He(5 3 D) + He + He* + e 

(119) He(5 3 F) + He -> He* + e 

Type 11: X*+e + e->X^ + X+ e 

(120) He* + e + e -* He(2 1 S) + He(l 1 S) + e 

(121) He* + e + e + He(2 3 S) + He(l 1 S) + e 

Type 12: X^+e + e+X^X+e 

(122) He* + e + He -* He(2 1 S) + He + He 

(123) He* + e + He -> He(2 3 S) + He + He 


26 

26 

26 
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24,25 

24,25 


24,25 

24,25 
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Type 13: X* + X* + X + + X + e 


(124) 

He(2 1 S) + He(2 1 S) He + + He + e 

40,41 

(125) 

He(2 1 S) + He(2 3 S) + He + + He + e 

40,41 

(126) 

He(2 3 S) + He(2 3 S) -*■ He + + He + e 

27,41,50 


When nitrogen is present the following reactions 

are included: 


Type 4: X+eJX + +e + e 


(127) 

N 2 + e J + e + e 

39 


Type 8: X* + e X + X 


(128) 

N* + e % N + N 

51 


Type 14: X + + N 2 N* + X 


(129) 

He + + N 2 X N* + He 

52 

(130) 

He 2 + N 2 i N 2 + He + He 

53 


* + 



Type 15 : X + -* + X + e 


(131) 

He(2 1 S) + N 2 + + He + e 

54 

(132) 

He(2 3 S) + n 2 ■* + He + e 

54 
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APPENDIX B 


GAIN COEFFICIENT 


The possible lasing transitions in atomic helium involve electronic state 
transitions. The gain coefficient for transition from upper level m to lower 
level n is 


Y(V) 


c 2 g. 


m 


8tt V 


2 n 

mil 


where G(v) is the shape factor. N and N 

r m n 

n, g and g are their degeneracies, A 
taneous emission from level m to level n, 
tion. 


N N 

[~ - G(V) (Bl) 

s m ®n 

the number densities of states m and 
is the Einstein coefficient of spon- 
and V is the frequency of the transi- 


The shape factor is affected by two processes, Doppler and pressure broad- 
ening. In general, for pressures greater than 30 Torr, pressure broadening is 
dominant. Thus, Doppler broadening is neglected in this analysis. For pres- 
sure broadening the shape factor G(v) is given by 


where 


G(V) = 


IT E V 

t St 


(B2) 


T T 1/2 

v . = 4 [2k c-- + -r)] z n 
st j m m. st t 

s t 


(B3) 


In equation (B3) , s represents the lasing gas and t represents other gases pre- 
sent in the system. When summing on t in equation (B2) , t can be equal to s. 

The quantities Z can be calculated from the Lennard-Jones potential 
st 


parameters of the colliding molecules 

Z 


60-62 


= d 2 fl< 2 ’ 2 V,) 

st st st 


(B4) 
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d = /d d~ 
st st 


* , * 

T = T/e 
st st 


* 

E st 


/* * fit * 

2 /e e /l I 

s t s t 

* * 

< T s + V 


2) 63 

The quantity ft'' * ^ (T) is obtained from the empirical curve fit° 

ft ( 2 , 2 ) (T) = + o > 5 2487 exp [-0.7732 T] + 

2.16178 exp [-2.43787 T] - 

, ,or „-4 m 0. 14874 . ,18.0323 

6.435 x 10 T sm (“q - 7533 " “ 7.27371). 


The required parameters for helium are 


60,64 


d = 2.576 x 10 ~ 8 cm 
s 


E = 10.22 °K 
s 


I = 24.586 ev . 
s 


(B5) 

(B6) 


(B7) 


(B 8 ) 
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Figure 1. Slab approximation of cylindrical geometry. 
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Figure 7. Helium excited states at 100 Torr, 2^P model. 
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Figure 8. Helium excited states at 760 Torr, 2 P model. 
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Figure 10. Helium excited states at 760 Torr, TM model. 
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Figure 11. Effect of excitation transfer reactions on 
excited states. 
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Figure 12 . F.ffect of associative ionization and excitation 
transfer rates for n - 4 on excited states. 


49 



3 3 S# 


• 3 S 


3 5 P».,3, 


# 3'p 


• 4 D 


4S* #4 3 d 
4 3 S» 


# 4 3 P 


<£=3.8x10 NEUTRONS/CM SEC 

P = 100 TORR 
N 2 = 0 PPM 
— E = 10 V/CM 


I I I I 


I I 


£ 


NASA-Langley, 1977 CR-2776 


